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We use computer simulations to study the relaxation dynamics of a model for oil-in-water mi¬ 
croemulsion droplets linked with telechelic polymers. This system exhibits both gel and glass phases 
and we show that the competition between these two arrest mechanisms can result in a complex, 
three-step decay of the time correlation functions, controlled by two different localization length- 
scales. For certain combinations of the parameters, this competition gives rise to an anomalous 
logarithmic decay of the correlation functions and a subdiffusive particle motion, which can be un¬ 
derstood as a simple crossover effect between the two relaxation processes. We establish a simple 
criterion for this logarithmic decay to be observed. We also find a further logarithmically slow relax¬ 
ation related to the relaxation of floppy clusters of particles in a crowded environment, in agreement 
with recent findings in other models for dense chemical gels. Finally, we characterize how the com¬ 
petition of gel and glass arrest mechanisms affects the dynamical heterogeneities and show that for 
certain combination of parameters these heterogeneities can be unusually large. By measuring the 
four-point dynamical susceptibility, we probe the cooperativity of the motion and find that with 
increasing coupling this cooperativity shows a maximum before it decreases again, indicating the 
change in the nature of the relaxation dynamics. Our results suggest that compressing gels to large 
densities produces novel arrested phases that have a new and complex dynamics. 

PACS numbers: 


I. INTRODUCTION 


In nature and in our daily life, many soft materials are 
formed due to the dynamical arrest of the constituent 
particles Usually they are labelled as gels if the 

particle density is low and as glasses if the density is 
large. However, the difference between these two states is 
at present not very well understood and therefore it is not 
always easy to distinguish them. Despite this difficulty 
quite a few features in this glass-gel cross-cover regime 
have been studied extensively. 

In dense glass forming liquids the slowing down of 
dynamics is related to the mutual steric-hindrance in 
the motion of the constituent particles. The dynam¬ 
ical properties of these glass formers are characterized 
by a stretched-exponential shape of relaxation functions 
[lj, or similarly by the anomalous, exponential tails in 
the van Hove distributions of particle displacements 3- 
Many of these dynamical features are described well by 
mode-coupling theory [lj. At even lower temperatures, 
the relaxation dynamics can be understood by means of 
the random first order transition theory Q. In these 
glass formers, the structure is usually close-packed for 
hard sphere or van der Waals type interactions and is 
accompanied by a super-Arrhenius increase of viscosity. 
Or, if there are covalent bondings, they form network-like 
structures and exhibit an Arrhenius increase in viscosity. 

On the other hand, chemical gels or rubbers are soft 


solids having random network structures @,i,0. The 
cross linking of the permanent bonds between the con¬ 
stituent monomers happens during the synthesis process 
inducing a vulcanization transition once the density of 
the links exceeds the percolation threshold js[| . Different 
static and dynamical properties in the vicinity of this 
transition have been studied, both using simulations and 
theoretical models (e.g., see 1 , 0 ). Physical gels are 
on the other hand low-density network structures with 
bonds that can be broken/re alig ned by thermal fluctua¬ 
tions within finite timescales [7|. One possible nonequi¬ 
librium path to physical gelation is via a thermal quench 
across the liquid-gas spinodal leading to dynamically ar¬ 
rested states HI- 141, that show complex aging phenom¬ 
ena 00- In general, these paths lead to spatially het¬ 
erogeneous structures. However, in recent times, consid¬ 
erable effort has been made to devise ways by which spa¬ 
tially homogeneous physical gels can be formed 0 , 0 - 
[0 . For such gel-forming systems, a wide variety of relax¬ 
ation functions have been reported: logarithmic 0-0], 
stretched 0 or compressed exponentials 0,0 . While 


theoretical models [26|, I32fl34j have been proposed to ac¬ 
count for such dynamical properties, they are at this time 
certainly not yet comprehensive. 

Of particular interest is the interplay between these dif¬ 
ferent arrest mechanisms, viz. gel and glass, since their 
competition can be used to engineer materials with novel 
functionalities 0,0-0] . Similar studies have been car¬ 
ried out in systems with competing lcngthscales 
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or interactions |H SU El ■ Yet, only few models do al¬ 
low to study the low density gel phase and high den¬ 
sity glass at the same time. Accessing this regime is, 
however, necessary for investigating the structure and 
dynamics at those intermediate densities where the gel 
transforms into a glass and vice versa. Here, we study a 
simple model with direct experimental relevance (l 8 | , and 
which permits us to traverse the density regime of inter¬ 
est and hence to study the interplay of different processes 
which lead to either gelation or glassiness. Furthermore, 
our model also allows us to tune the lifetime of bonds, a 
feature that is usually not present in other models (for 
example, see the recent work [42|] ). On one hand, this 
facilitates a wider exploration of the relaxation dynam¬ 
ics of such model physical gels, but also allows to disen¬ 
tangle the origin of the apparently anomalous relaxation 
observed in these systems. In fact, mode-coupling theory 
predicts, e.g., logarithmic relaxation whenever two differ¬ 
ent arrest lines meet, one gel-like and another glass-like, 
and relates it to an underlying higher-order singularity 
in the theory I33j|. The versatility of our model, and 
in particular the possibility of tuning at will the bonds 
lifetime, allows us to explore the different mechanisms 
at play and hence to elucidate the interplay of various 
lengths and time-scales [ 20 ]. 

The paper is structured as follows. In Section II we 
explain the details of our model transient-network fluid, 
together with the numerical schemes used to simulate 
its dynamics. The phase diagram of the model system 
and its structural properties are discussed in Sections III 
and IV, respectively. In Section V, we analyze in de¬ 
tails the fluid’s dynamics, quantified by the mean squared 
displacement and the incoherent scattering function. For 
that, we follow different routes across the phase diagram, 
which allows us to clearly understand the interplay be¬ 
tween the gel and the glass regimes. Finally, the dynami¬ 
cal heterogeneities which characterize the slow dynamics 
in the gel and glass phases are studied in Section VI, 
followed by a summary of our results and a broader per¬ 
spective in Section VII. 


II. MODEL AND DETAILS OF SIMULATION 


Our model system is a coarse-grained representa¬ 
tion G2 of a transient gel which has been studied in 
experiments [18|. In this system, an equilibrium low- 
density gel is obtained by adding telcchclic polymers to 
an oil-in-water microemulsion. Since the polymer end- 
groups are hydrophobic, the polymers effectively act as 
(attractive) bridges between the oil droplets they con¬ 
nect. The strength, lengthscale and typical lifetime of 
these bridging polymers can be controlled at will. Denot¬ 
ing by Cij the number of polymers connecting droplets 
i and j, we have established in Refs. [Hill , that the 
following interaction is a reasonable coarse-grained rep¬ 
resentation of this ternary system: 


V = ei 


j>i 
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(1) 

The first term is a soft repulsion acting between bare 
oil droplets, where a-ij = (cq + Uj)/ 2 , eq is the diam¬ 
eter of droplet i, and is the distance between the 
droplet centers. The second term describes the entropic 
attraction induced by the telechelic polymers, which has 
the standard “FENE” (finitely-extensible nonlinear elas¬ 
tic) form known from polymer physics 0 , Vfene(Uj) = 
ln(l — {rij — <Jij) 2 /£ 2 ), and accounts for the maximal 
extension £ of the polymers. The last term introduces 
the energy penalty eo for polymers that have both end- 
groups in the same droplet. The most drastic approxima¬ 
tion of the model m is the description of the polymers 
as effective bonds between the droplets, which is justi¬ 
fied whenever the typical lifetime of the bonds is much 
larger than the timescale for polymer dynamics in the 
solvent [18]. Thus, for exploring the different properties 
of such a model, the relevant variables are the droplet 
volume fraction (f> and the number of polymer heads per 
droplet R E3- In order to describe the dynamics of the 
system, we use a combination of molecular dynamics to 
propagate the droplets with the interaction ©, and lo¬ 
cal Monte Carlo moves with Metropolis acceptance rates 
t - iTriki 11 i n [ 1 , exp(A V/ksT)) to update the polymer connec¬ 
tivity matrix Cij , where AV is the difference in poten¬ 
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of the system for the two bond configura- 
19j . Thus Tii n k is the timescale governing the 
renewal of the polymer network topology. In order to 
prevent crystallization at high volume fractions (which 
would be the case for the monodisperse model studied 
earlier 0 ), we use a polydisperse emulsion with a flat 
distribution of particle sizes in the range cq S [0.75,1.25] 
(having a mean diameter a = 1). The units of length, en¬ 
ergy and time are respectively a, and m/ei where 
m is the mass of the particles. The space of control pa¬ 
rameters is quite large. Therefore we set £ = 3.5 a as 
measured in experiments jl 8 }, T = 1, and eo = 1 and 
62 = 50, and vary the remaining parameters {<fi, R , ru n k}. 
These choices for the parameter values leads to a phase 
diagram which is similar to the one obtained in experi¬ 
ments. 

Our numerical simulations are done for a three dimen¬ 
sional system of N = 1000 particles. The equations of 
motions of these particles are integrated using a velocity 
Verlet scheme with a time step of St = 0.005. Here, most 
of the results are reported for ru n k = 10 2 , although we 
also explore other values of ri; n k: 1 , 10 , 10 3 , 10 4 to illus¬ 
trate some of the dynamical features of the system. At 
each volume fraction <(>, we first equilibrate the system of 
particles without any links (R = 0). Once equilibrated, 
bonds are introduced corresponding to the required value 
of R and then the system is again equilibrated to obtain 
the proper distribution of bonds per particle. Since the 
structure of the network is independent of the choice of 
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FIG. 1: (Color online) Phase diagram of the system obtained 
by varying R and <j). The diamonds correspond to the coex¬ 
istence region between gas and liquid, the circles to the sol 
phase and the squares to the gel phase. The horizontal dashed 
line corresponds to onset of percolation, the dot-dashed line 
indicates the phase-coexistence boundary and the thick line 
marks the predicted glass line. Snapshots of typical configu¬ 
rations in each phase have been published elsewhere, see Fig. 
3 in 0. 

Tink El , we use a small value of rii n k = 1 to expedite the 
equilibration process. Subsequently, data is generated 
by continuing the simulations with different Tii n k values 
when required and the averages are typically calculated 
over 100 different time origins. We also do simulations 
for the case when the bonds between two particles are 
completely frozen. In order to do a proper sampling of 
the network configurations for this situation, we use 6 
initial configurations (positions, connectivities) from the 
simulations with a finite rii n k as initial inputs for subse¬ 
quent evolution of the particle positions using molecular 
dynamics with the connections now permanently fixed. 


III. PHASE DIAGRAM 

We begin by summarizing our earlier findings for the 
phase diagram (shown in Fig. [T]) for this model. If the 
number of polymers is small (i.e R < 2), the system is in 
a simple liquid phase (the sol) at small values of (f). In this 
regime, the distribution of connectivities per particle is 
just an exponential |l9j. With increasing (j), the dynam¬ 
ics in the sol regime becomes slow and one eventually 
enters a glassy phase at large (j), characterized by a very 
strong increase of the timescales for structural relaxation. 
If the number of bonds is large and (f> is small, phase sep¬ 
aration is observed due to strong attractions between the 
droplets. Here, the distribution of connectivities becomes 
bimodal, with one of the two peaks corresponding to a 
well-connected liquid and the other to free particles. For 
intermediate values of R the system is in a gel phase. In 
this region of the phase diagram the particles are con¬ 


nected together and form a percolating cluster and the 
spatial density is homogeneous. Here, the connectivity 
distribution is peaked around the average value and has 
an exponential tail. Note that this gel is an equilibrium 
phase since the polymer network is constantly rewired 
on the timescale Ti; n k. However, if we go to large enough 
volume fractions, we observe again a glassy system for 
all connectivities, with the corresponding divergence of 
relaxation timescales. 

At low <f >, in the gel phase, the main slow relaxation 
process is related to the connections between the droplets 
by means of the polymers and the timescale associated 
with its reconfiguration. At low R , as the system becomes 
glassy at large (j), the origin of the slow dynamics is the 
steric hindrance caused by the caging of each particle by 
its neighbors. In the region where both R and 4> act as a 
source for slow dynamics, we have shown that the generic 
relaxation process has three steps but with proper tun¬ 
ing of the two relaxation timescales, one can also obtain 
logarithmic decays of the relaxation function [2Cj . 

In the following we will discuss in details the interplay 
between these two relaxation processes and show the con¬ 
sequences on the nature of different dynamical quantities 
as we move around in the phase diagram. 


IV. STRUCTURE 



FIG. 2: (Color online) Main panel: Structure factor S(q) 
computed for the particles at a volume fraction of <f> = 0.70 for 
different connectivities R = 0,4, 8, and 12. Inset: S(q) at (f> = 
0.50 for R = 0,12, 30, and 50. The dashed line corresponds 
to q~ 4 . 

Before we discus the different dynamical properties of 
the system, we briefly look at its structure. In Fig. [2] we 
plot the static structure factor, S(q), for a system that is 
dense, (j) = 0.70, varying the connectivity R. The general 
shape of S(q) is very similar to the one of a simple liquid 
and hence we can conclude that the system is homoge¬ 
neous. Also, we see no significant dependence of S(q) on 
R, thus showing that at high densities the structure is 
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mainly governed by steric hindrance. For intermediate 
densities, however, we do note the weak dependence on 
R in the regime of small wave-vectors, illustrated in the 
inset of Fig. [2] where we show S(q) for <j> = 0.50. This 
is at a volume fraction at which for large R the system 
approaches the co-existence region, and hence one starts 
to see the emergence of a power-law behavior at small 
values of q with increasing i?; the data for R = 50 can 
be approximated by q -4 (which is expected in proximity 
to phase co-existence [43}). 


V. RELAXATION DYNAMICS 
A. Dependence on volume fraction <j> 

We now characterize the dynamical proper¬ 
ties of the system by focusing on two quanti¬ 
ties: (i) the mean squared displacement, defined 

as A 2 {t) = (|ri(t) - r,(0)| 2 ) and (ii) the 

self-intermediate scattering function, defined as 
F s {q, t) = (jf E j ex P (*q ' - ij-(0)])>. Here r t (t) is 

the position of particle i at time t, q is the wave-vector, 
and (.) corresponds to the ensemble average. 

For increasing volume fraction </> = 0.50,0.61,0.70, 
we discuss simultaneously the data for A 2 (t), shown in 
Fig.[3j and F s (q, t), shown in Fig. [31 computed at a wave- 
vector value q = 6. Thus, the measured F s (q,t) probes 
the relaxation dynamics on length scales that are slightly 
larger than the average particle diameter (the peak in the 
structure factor S(q) occurs at q m 7.3). 

We start in the pure gel phase (cj> = 0.50). The mean 
squared displacement, Fig.[3jr, shows that the increasing 
number of bonds restricts the motion of the particles in 
that for R > 0 we see the emergence of an intermedi¬ 
ate regime which develops into a well-defined plateau at 
R = 10. The height of this plateau depends significantly 
on R , showing that this cage motion is directly related to 
the transient bonds between the particles. The self inter¬ 
mediate scattering function, Fig. 0Ji, shows for small R a 
very rapid decay. This changes in that for R around 4-6 
a plateau develops at intermediate and long times, the 
height of which depends strongly on R. The presence of 
this increasing plateau height, which is reminiscent to the 
so called type-A transition of mode-coupling theory , 
indicates that the relaxation mechanism is changing: For 
small R the motion of the particles is only weakly slowed 
down by the presence of the bonds, which typically break 
on the time scale of m n k. However, for larger R breaking 
a few bonds is not enough to allow the particles to move 
since the remaining bond still allow to maintain the par¬ 
ticle inside its cage. Hence this makes that at large R the 
relaxation dynamics does not depend very strongly on R 
anymore. This effect is seen in Fig. [5] where we show 
the diffusion constant of the particles, D, (as obtained 
from the mean squared displacement at long times) as 
a function of R. For small R, D shows a rather strong 
R— dependence, whereas for R > 5 this dependence be- 



FIG. 3: (Color online) Variation of mean squared displace¬ 
ments A 2 (t) with changing R for (a) cf> = 0.50, (b) tj> = 0.61, 
and (c) <j> = 0.70, using a bond lifetime of rn n k = 10 2 . 


comes weaker. Below we will discuss the R —dependence 
of the relaxation time in more detail. 

Next, we look at the data for an intermediate density, 
viz. <j) = 0.61, and in Fig. [3)o we show the correspond¬ 
ing A 2 (f). Like for (f> = 0.50, the longtime diffusion de¬ 
creases with increasing R, and the R —dependence of the 
diffusion constant shows again a break at around R « 5 
(see Fig. 0. For this value of <f> we observe, however, 
for R > 2, at intermediate times a shoulder in A 2 (i). 
This shoulder, clearly visible for R = 4, is related to 
the presence of the bonds that lead to a caging of the 
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FIG. 4: (Color online) Variation of the self intermediate scat¬ 
tering function with changing R. (a) <j> = 0.50, (b) <j> = 0.61, 
and (c) <j> = 0.70, using ru n k = 10 2 . Note the change of time 
span in different panels, from a maximum time of 10 4 in (a) to 
10 6 in (c). The inset in panel (c) shows the collapse of corre¬ 
lation functions for 0 = 0.70 after rescaling by the relaxation 
time, similar to time-temperature superposition principle. 


particles on the length scale related to £, the maximum 
extension of a bond. Thus the hint of the short-time 
plateau (when A 2 {t) ss 0.1) and again one at later time 
(when A 2 (f) w 1) reflects the presence of the two dif¬ 
ferent mechanisms for constraining particle motion, viz. 
local steric hindrance and the network bonds. Since each 
type of caging leads to a plateau in the intermediate scat- 


Q 


FIG. 5: (Color online) Variation of diffusion constant D with 
R, for different <j> (shown in Fig. [3]). The dotted line, corre¬ 
sponding to an exponential function, is drawn as a guide to 
the eye. 


tering function (2C|, the existence of the two competing 
mechanisms makes that at intermediate times F s (q,t ), 
shown in Fig. [Dy has a very slow, almost logarithmic, 
decay, if R is around 6. More details on this particu¬ 
lar case are give in the context of Fig. [FT] If one com¬ 
pares the data for R = 10 at the two volume fractions 
(f> = 0.50 and <f> = 0.61, we recognize that the height of 
the plateaus in F s (q,t) increases with <f> from which one 
can conclude that the proximity of the particles leads to 
increased tightening of the cage. 

If the density is increased further to (f = 0.70, the re¬ 
laxation dynamics becomes strongly dominated by the 
steric hindrance mechanism. Already for R = 0 one sees 
a weak plateau in A 2 (i) and its length grows rapidly 
with increasing R without changing much its height (see 
FigEt)- This is the typical behavior of simple glass¬ 
forming liquids [44j . At the same time the self inter¬ 
mediate scattering function shows the growth of a shoul¬ 
der with finite height and this height depends again only 
weakly on R, Fig. 0J:. In contrast to the case at lower 
densities, here the shape of the correlator is basically in¬ 
dependent of R. This is demonstrated in the inset of 
Fig. where we plot F s (q,t) as a function of i/r, with 
the relaxation time r defined by F s (q,r) = 0.03. The 
fact that this presentation of the curves leads to a nice 
master curve shows that we have for this system a time-i? 
superposition, in analogy to the time-temperature super¬ 
position found in simple glass-forming systems [l] . 

Despite the qualitative changes seen in A 2 (t) and 
F„(q,t) if <j> is increased, the I?—dependence of the diffu¬ 
sion constant for (f> = 0.70 is very similar to the one seen 
at lower densities, see Fig. [5j Also for this high value 
of </> we see that this dependence is relatively strong at 
small R and becomes weaker if R > 6. Hence the fact 
that at low R just few bonds have to be broken in order 
to allow a particle to move whereas at high R this is not 
a sufficient condition, is reflected in the J?—dependence 
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FIG. 6: (Color online) F s (q,t) for R = 6.3 ,<f> = 0.61, a) 
Tiink = 10 3 b) ri ink = 10 2 and c) ri ink = 10. 

of D at all <f>. 

B. Varying the bond life-time Tiink 

We will now explore further the interplay between the 
two processes leading to slow relaxation, i.e. the nearest- 
neighbor caging and the constrained motions due to the 
polymer bonds. Since the relative importance of these 
two processes depends on the lifetime of the polymer 
bonds, we will in the following vary this lifetime and con¬ 
sider values of rn n k = 10,10 2 , and 10 3 , and at fixed state- 


point of R = 6.3, ^ = 0.61, where logarithmic decay in 
the time-correlation function is observed. We study how 
the shape of F s (q,t) changes with varying Tn n k and will 
relate this to the interplay between the two processes. 
This is done for different values of wave-vector in order 
to see the how relaxation timescales vary over different 
lengthscales. 

We begin by looking at the case of the large Ti in k = 10 3 , 
i.e. when the polymer bonds hinder the motion on a 
time scale longer than the steric hindrance effect, see 
Fig. [HJi. The correlation function reflects three differ¬ 
ent relaxation processes which can be seen for all value 
of q. Initially, the particles rattle inside the cage, result¬ 
ing in partial relaxation of the correlation function on a 
timescale Tp ps 1. Later on, the particles escapes from the 
cage of neighboring particles, (which for this value of <f> 
not very pronounced), but the relaxation process is then 
held up by the polymer bonds. Eventually, the polymer 
network rewires on a timescale which is proportional to 
Tiink = 10 3 and the particles start the final relaxation pro¬ 
cess. The height of the plateau in F s (q,t) increases with 
decreasing q, which is the typical behavior for a glassy 
system [44], E]. However, for the range of wave-vectors 
explored, the final timescale for decay of F s (q,t) depends 
only weakly on q. Thus, over these length scales, the re¬ 
laxation process is determined by the reconfiguration of 
the network. However, at larger length scales, one can ex¬ 
pect that hydrodynamic effects will eventually dominate 
and this will then determine the relaxation timescales. 

For intermediate values of the bond lifetime, rii nk = 
10 2 in Fig. [6 Jd, the final decay has moved to shorter times 
and makes that now the interplay between the two pro¬ 
cesses results in a logarithmic decay of F s (q,t) as dis¬ 
cussed elsewhere (20|. This logarithmic dependence is 
seen for a range of q values (see Fig. [6jo), with the time- 
window over which it exists decreasing with decreasing 
q. This implies that this form for the correlation func¬ 
tion occurs only for specific combination of the relaxation 
timescales of the two processes of steric hindrance and 
eventual network relaxation. At sufficiently small q , hy¬ 
drodynamics makes that the relaxation becomes so slow 
that the second relaxation step is no longer visible out 
and hence the logarithmic t —dependence is no longer ob¬ 
served. Finally we mention that the logarithmic shape in 
the relaxation function as discussed here is not related to 
any underlying higher order mode-coupling transition, in 
contrast to the case of certain colloidal systems for which 
similar relaxation functions have been observed [27j . 

If we set rii n k to a small value, this three-step relax¬ 
ation can no longer be observed, as is seen for the case 
of Tii n k = 10 in the bottom panel of Fig. [6] For all q’s 
the curve show a (seemingly) simple two step relaxation, 
since the third step (related to the bonds) starts already 
when the second step (related to steric hindrance) is not 
yet completed and hence the two processes become com¬ 
pletely mixed in time. Below we will briefly come back 
to this effect. 
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C. When the bonds are permanent 

A useful way to check the influence of the polymer 
bonds on the relaxation dynamics is by comparing the 
relaxation functions for the case when the bonds are per¬ 
manent, i.e. 7ii n k = oo, to those when the lifetime is finite. 
In Fig. 0 we do this comparison for different connectivi¬ 
ties (R = 4,6, and 8) and different values of </>. 

In Fig. [7H, we show F s (q = 6, t) for the case of R = 4 
with the two different bond lifetimes Tii n k = 10 2 and 
Tiink = oo. For cj) = 0.50, 0.60, and 0.65, the timescales 
for overcoming the steric hindrance are the same for both 
lifetimes. However, for all </>, we find that the correlation 
function for Tii n k = oo shows a plateau at long times 
(not visible in this plot), which is due to the fact that 
the frozen bonds in the percolating gel-network prevent 
the complete relaxation of the system. In contrast the 
curves for the finite rn n k vanish at long times. For small 
and intermediate (f> the two sets of curves are very sim¬ 
ilar, indicating that the presence of a few bonds does 
not change the dynamics significantly. Only for tf> = 0.7 
one sees a substantial difference in that the correlator 
for the permanent links decays slower than the one with 
Tiink = 10 2 . It is reasonable that these differences are 
noticeable for times somewhat longer than 10 2 , i.e. the 
time scale of rii n k- 

If we increase the connectivity to R = 6 (Fig. EH>), we 
see that the behavior is qualitatively similar to R = 4 in 
that for all values of 4> the two curves track each other 
up to times around 10 2 , i.e. the time of the finite rii n k- 
For larger times the correlators for Tn n k = 10 2 decay to 
zero whereas the ones for rii n k = oo show at long times 
a marked plateau. The height of this plateau depends 
now more strongly on <f> than it was the case for R = 4, 
showing that if R is increased the life time of the bonds 
becomes more influential. This is reasonable since it is 
related to the general observation that in glass-forming 
system small changes influence the relaxation dynamics 
increasingly more the slower the dynamics is. We also 
note that for <f> = 0.70 the correlator for the permanent 
bonds becomes very stretched. This sluggish relaxation 
might be related to the fact that for this value of R there 
are, in addition to the percolating cluster, clusters of dif¬ 
ferent sizes (see Ref. [H] for typical distributions), thus 
giving rise to relaxation dynamics that spans many or¬ 
ders of magnitude in time and hence to a very stretched 
average correlation function. The stretching of the cor¬ 
relator for the frozen bonds could, however, also be due 
to the fact that these different clusters hinder each other 
resulting in the overall slowdown of the dynamics [46[ . 

Next, we increase the number of bonds even more, viz. 
R = 8, as shown in Fig. [7J:. We see that for the case 
of permanent bonds the height of the asymptotic plateau 
has increased strongly in comparison to the case of R = 6. 
As a result the correlator for </> = 0.70 shows only a neg¬ 
ligible decay of the correlation function if the bonds are 
permanent. The motion is so constrained by these bonds 
that the height of the asymptotic plateau, caused by the 



FIG. 7: (Color online) F 3 (q,t) for different cf> and R = 4,6, 
and 8 when the bond lifetime is rii n k = 100 (lines with sym¬ 
bols) and when rij n k = oo (lines). The wave-vector is q = 6.0. 


permanent bonds, becomes comparable to the one re¬ 
lated to the steric hindrance. Also for tfi = 0.5 the height 
of the second plateau has increased so much that the 
relaxation from caging is now completely masked. How¬ 
ever, for the intermediate values of <j), one does notice a 
difference between the two plateau heights and the corre¬ 
lations functions decay in a very stretched fashion from 
one to the other. In fact, the decay is so slow that the 
time-dependence is seen to be logarithmic (nearly for five 
decades in the case of 4> = 0.65). Note that this loga- 
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FIG. 8: (Color online) (a),(b): Variation of F a (q, t) with wave- 
vector q for R — 2 and 4 at <f> = 0.60 with rii n k = oo. (c): 
Height of plateau, fq N , of F„ (q. t) at long times as a function 
of q for R = 2, 4, and 6. The dashed lines are fits of the form 
exp (— q/Q with the corresponding £ = 0.55, 0.99, and 2.27. 
(d),(e): vs. q for different values of <j> for R = 4 and 6. 


rithmic decay is due to the heterogeneous relaxation of 
the floppy clusters of frozen bonds, which we will discuss 
later in further detail. Thus this mechanism is different 
from the one leading to the logarithmic relaxation seen 
in Fig. |4] which was due to a unique combination of the 
two finite relaxation timescales. 

We now study the Soppiness of this network of particles 
connected by the permanent bonds by probing the wave- 
vector dependence of the relaxation functions F s (q,t). 
In Figs. [5Ji and b, we show for (j) = 0.60 the variation 
of F s (q,t) for R = 2 and 4, i.e. in the region of the 
phase diagram where gelation sets in. The height of the 


plateau at long times, also called non-ergodicity param¬ 
eter, is a measure for the stiffness of the network on the 
length scale q considered. Comparing the two panels we 
recognize that, for a given q, the height of the plateau in¬ 
creases with increasing R. Denoting this height by /™, 
we show in Fig. [5J: that /™ shows basically an expo¬ 
nential decrease in q with a slope that decreases rapidly 
with increasing R. That /™ decreases with increasing 
q is of course reasonable since on small length scales the 
particles have more leeway to flop around than on large 
length scales. Note, however, that this exponential de¬ 
pendence is in contrast to the one found for the height of 
the plateau due to the steric hindrance, the latter being 
basically a gaussian function [45| . Since in the represen¬ 
tation of Fig.jSJ: such a gaussian dependence is given by a 
parabola, we see that such a curve will intersect the one 
for at a certain value of q x . For q < q x the plateau 
due to the steric hindrance is above /^ N , thus making 
that one observes two plateaus in the correlator. How¬ 
ever, for q > q x the plateau at long times is higher than 
the steric one, thus making that the latter one will be 
completely masked by the former and thus the correlator 
will show only one plateau. 

We have also studied how the q —dependence of /^ N 
changes with the volume fraction and in Figs. [SJI and e, 
we show for R = 4 and 6, respectively. For both 
cases we see that the rigidity of the network at large 
scales, i.e. small q , is not affected by the volume fraction 
which is not surprising. For R = 4 we see that this is also 
true at small length scales whereas for R = 6 we note a 
significant (^—dependence if q is large. This difference is 
likely related to the fact that the system with R = 6 is 
much more sluggish than the one for R = 4, see FigJT] 
and hence small changes (here in (j>) will have a stronger 
impact on the dynamics. 

Finally we disentangle the dynamics of the particles 
that belong to the percolating cluster from the one that 
are not attached to it. In the following discussion these 
are referred as clustered and non-clustered particles, re¬ 
spectively. The objective is to clarify the respective con¬ 
tributions to the different dynamical quantities that we 
have discussed above. We do this comparison for an in¬ 
creasing number of connections I? at a fixed (large) vol¬ 
ume fraction of (f> = 0.65. In Fig. [9] we show the data for 
A 2 (i) and F s (q,t) for these two families of particles. 

The mean squared displacement of the clustered par¬ 
ticles shows after the ballistic regime at short times a 
shoulder that is related to the cage of the steric hindrance 
(Fig. Et 1 ) - This localization is, however, only temporary 
and is followed by a further increase of A 2 (t). Only at 
longer times A 2 ( t ) saturates at a height that depends on 
R. We note that the approach to this asymptotic height 
becomes increasingly slow with increasing R and in fact 
for R = 8 the time dependence is close to logarithmic 
and does not end within the time window of our simula¬ 
tion. This behavior is also seen in the self intermediate 
scattering function F s (q,t) (Fig.^r). At short times the 
correlator decays quickly onto a plateau (not very pro- 
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FIG. 9: (Color online) Relaxation dynamics of the clustered 
(dashed lines with open symbols) and non-cluster particles 
(full lines with filled symbols) for the case of frozen bonds (i.e 
Tank = oo), at R = 4, 6, and 8 for volume fraction <j> = 0.65. 
(a): Mean squared displacement A 2 (t). (b) Self intermediate 
scattering function F a (q,t) for q = 6. 

nounced) before the relaxation of the steric hindrance 
starts. For R = 4 and 6, this process ends in that the 
correlator reaches the final plateau (which is given by 
/^ N discussed above). However, for R = 8 the final de¬ 
cay is so slow, again compatible with a logarithmic time 
dependence, that we do not see the asymptotic behavior. 

Finally we look at the motion of the non-clustered par¬ 
ticles and compare it with the one for R = 0. From 
Fig. [9ji we recognize that at this volume fraction also 
these particles are slowed down by the cage effect in that 
one sees for all values of R a shoulder in A 2 (t) at time 
t ~ 1. For R = 0 the A 2 (t) shows then immediately 
the diffusive behavior, i.e. it is proportional to t. How¬ 
ever, if R is increased, the t —dependence of A 2 (t) for the 
un-clustered particles follows first the one of the clus¬ 
tered particles. Only once the latter starts to reach the 
plateau discussed above, do the former cross over to the 
diffusive behavior. Hence we can conclude that before 
this crossover the relaxation dynamics of the two pop¬ 
ulation of particles are strongly coupled. This result is 
reasonable because in order to move, the un-clustered 
particles have to explore the holes within the percolating 


cluster formed by the particles which are permanently 
linked. 

Also the self intermediate scattering function of the 
non-clustered particles tracks the one of the clustered 
particles at short times (Fig. [9jo). However, once the 
latter starts to show at long times a plateau that has a 
significant height, the two correlators differ strongly since 
the one for the non-clustered particles decays to zero at 
long times. From the graph we also see that for R = 8 the 
correlator is extremely stretched and shows almost a log¬ 
arithmic t— dependence. This very slow decay indicates 
that the mobile clusters can move around the percolat¬ 
ing cluster only with great difficulty, a behavior that is 
similar to the relaxation dynamics of particles moving in 
random porous media [13, SH . It is also interesting that 
for the highest R the mean squared displacement shows 
for the last two decades in time a nice diffusive behavior, 
whereas F s (q, t ) is far from having decayed to zero. This 
apparent contradiction is related to the fact that A 2 (t) 
is dominated by the particles that move relatively fast 
(i.e. they are in the small clusters) whereas F s (q,t) is 
dominated by the slowly moving particles (i.e. the large 
clusters). This cluster-size dependent dynamics leads to 
a so-called dynamical heterogeneity and in Section fVTl we 
will discuss this phenomenon in more detail. 

Although we show in Fig. 0 the comparison between 
the dynamics of clustered and non-clustered particles for 
the case that Tii n k is infinitely large, it is evident that 
for a very large but finite value of rii n k the relaxation 
dynamics will be very similar. Hence if, e.g., Tn n k is on 
the order of 10 5 , basically none of the shown curves will 
change significantly and thus the conclusions drawn from 
Fig. [5] will be apply also for such value of Tn n k. 


D. Relaxation timescales 

We now investigate how the two different relaxation 
timescales, one due to breaking of local cages and the 
other due to the reconfiguration of the network-bonds, 
vary with the volume fraction and the connectivity R. 

To start we consider the case of structural relaxation 
related to the steric hindrance. In order to avoid that this 
relaxation process is influenced by the one of the network, 
we consider the case in which the latter is completely 
suppressed, which can be achieved by choosing rii„k = 
oo. In the following we will study the relaxation times 
associated with the intermediate scattering function for 
wave-vector q = 6. As discussed above, this correlator 
shows at long times an asymptotic plateau the height 
of which, /™, depends on R and (f >. To take this into 
account we define the relaxation time tsh, to be the time 
at which F s (q, t ) — J™ = 0.03. The evolution of tsh with 
</> is shown in Fig. [TOh . for different value of R. 

We see that in the absence of any bonds, i.e. for R = 0, 
we have the usual slowing down of dynamics with in¬ 
creasing (f>. The (^—dependence of the relaxation time 
can be fitted by a Vogel-Fulcher-Tammann-law of the 
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FIG. 10: (Color online) (a) Main panel: Plot of tsh vs. <j> for 
different values of R. Inset: Same data scaled by tsh (?,</> = 
0.5). (b) Main panel : tpn vs. R, for different values of tj). 
Inset: Same data as in the main panel but normalized by the 
relaxation time for R = 0, leading to a master curve. Note 
the double logarithmic scale in this plot. 

form tsh ~ exp [ A /(( f> c — with /3 ~ 1, from which 
we can estimate the volume fraction (f> c at which the re¬ 
laxation times would diverge. If we increase the number 
of bonds among the particles, we see that tsh increases. 
That this increase is not just a constant (I?—dependent) 
factor but depends also on (j> is demonstrated in the in¬ 
set of the figure where we have normalized the relax¬ 
ation times to its value at <j> = 0.5. Using the Vogel- 
Fulcher-Tammann-law we can thus extract from these 
data the I?—dependence of (f> c , which can be considered 
as a proxy for the I?—dependence of the glass transition 
temperature. This 4> C (R) line is included in Fig. Q] as 
well and we see that it has a weak negative slope with 
4> C (R = 0) = 0.847 and 4> C {R = 6) = 0.808. Thus we 
can conclude that the glass transition as induced by the 
steric hindrance mechanism does not depend strongly on 
the value of R. However, since the prefactor of the Vogel- 
Fulcher-law does depend strongly on R , we can conclude 
that a line of iso-relaxation time will bend significantly 
more. 

Next, we investigate how the timescale for the full re¬ 
laxation of the network of particles depends on the aver¬ 
age number of bonds, R. For this we define a relaxation 
time using the t— dependence of F s (q = 6, t) with a very 
short lifetime for the bonds (Tii n k = 1). tpn is then de- 



FIG. 11: (Color online) F s (q,t) vs. t for R = 6.3, (j> = 0.61, 
and different Ta n k. The inset shows A 2 (t) for the same param¬ 
eters. The dotted lines in the main figure mark the plateaus 
in F s (q,t ): /® H and /^ N , which respectively correspond to 
the localization lengths Isn and Ipn indicated in the inset. 


fined via F s (q = 6,tpn) = 0.03 and the R —dependence 
of this relaxation time is shown in Fig. [TUb . The figure 
shows that for small and intermediate values of R the /in¬ 
dependence is independent of the volume fraction in that 
the curves for the different values of (j> seem to be just 
shifted vertically. That this is indeed the case is demon¬ 
strated in the inset where we show that a plot of the 
same data, but now normalized to the relaxation time 
for R = 0, gives a master curve. The main panel shows 
that at small concentration of bonds the R —dependence 
is close to an exponential, a result that is likely related 
to the fact that increasing R leads to a tighter cage for 
the steric hindrance and hence a slower dynamics. How¬ 
ever, for intermediate values of R the curve tpn start to 
bend over towards a weaker R— dependence. In order to 
investigate this effect better we have carried out simula¬ 
tions for (f> = 0.50 at very high values of R: 16, 30, and 
50. We find that in this regime the relaxation time fol¬ 
lows closely an exponential (see main Fig. 110b ). but with 
an exponential scale that is smaller than the one seen at 
small R. Note that this dependence implies that there 
is no singularity in the relaxation dynamics at any finite 
value of R, at least for this volume fraction. Instead the 
dynamics shows a behavior that is similar to an Arrhe¬ 
nius law in that the barrier for the relaxation depends 
only on the number of bonds between the particles. This 
observation is in agreement with earlier simulations for 
equilibrium gels [5l| . 

We conclude this discussion by using the different 
relaxation timescales to develop a criterion that tells 
whether or not the correlator F s (q 1 1 ) will show the loga¬ 
rithmic t— dependence discussed in the context of Fig. 0b- 
For this we define /“ as the height of the plateau associ¬ 
ated to the a-process, with a £ {SH, PN}, and recall that 
Tp, tsh and tpn oc Ti; n k are the timescales for the three 
different relaxation processes described above, namely 
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the rattling inside the cage, the escape from the local 
cage, and the network renewal process. We can now de¬ 
fine a simple criterion for the anomalous logarithmic re¬ 
laxation to be observed by requiring that the slope of the 
two segments defined by pairs [(In rg, /® H ), (In rgn, /™)] 
and [(InrsH,/™), (hirpN,0)] is the same, see Fig. |TT] 
It is known that for glass-forming systems the plateau 
related to steric hindrance is a Gaussian functions of the 
wave-vector q, /® H ~ exp(—g 2 £g H ) 0. We showed ear¬ 
lier, see Fig. [Ht, that for the PN-process, /™ has an 
exponential shape at large q (the regime corresponding 
to anomalous logarithmic behavior). Putting these ele¬ 
ments together we thus obtain 


(qt pn - q 2 is h) = In 


1 ln(r SH r^*) 

ln(lp N 7sH ) 


( 2 ) 


This equation relates a purely structural observable 
which depends on the competition between two local¬ 
ization lengthscales with dynamical information encoded 
in the different relaxation times, predicting a precise 
connection between structure and dynamics whenever 
anomalous logarithmic relaxation is to be expected. 


VI. HETEROGENEITIES IN DYNAMICAL 
PROPERTIES 


Typical glass formers show a heterogeneous dynamics 
of the particles when the system is increasingly super¬ 
cooled. It reflects the broad distribution in the timescales 
for local structural relaxations (HI- On the other hand, 
we have reported earlier E3 that also the dynamics in 
the gel phase can be heterogeneous: The particles in the 
percolating cluster and those that are unattached have 
different nobilities till the timescales at which the bonds 
in the network are reconfigured. In the following we will 
thus explore how these two different sources of hetero¬ 
geneity interact as we increase the density of particles as 
well as the number of connectivities in the system. 

We use two different measures of dynamical hetero¬ 
geneity. The first one is the non-Gaussian parameter 
« 2 (t), defined as a 2 (i) = 3(r 4 (t))/5(r 2 (t)) 2 — 1, where 
(r 2 (t)) and (r 4 (t)) are the second and fourth moments of 
the distribution function of single particle displacements 
G s (r,t) = N- 1 J2i( s ( r ~ l r *:(i) “ d( 0)|)), i.e. of the self 
part of the van Hove function. A non-zero value of 0:2 (t) 
quantifies the extent of deviation from a Gaussian shape 
for G s (r,t) HH. Note that G s (r,t) is Gaussian at short 
times, i.e. in the ballistic regime, and again at long times 
when the particles are diffusive. 

The second measure is the dynamic susceptibil¬ 
ity t), computed via the fluctuations of time- 

correlation function: Xi(q,t) = N[(F^(q,t)) — (F s (q,t)} 2 ]. 
It is designed to capture the spatiotemporal correlations 
of particle nobilities and provides, from its peak value, an 
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FIG. 12: (Color online) Time dependence of Q 2 for different 
values of R (as marked in the graphs) for rn n k = 10 2 . (a): 
4> = 0.50 in double-logarithmic representation and (b): cj> = 
0.70. 


estimate of the dynamic correlations [5(1 [HI • Such func¬ 
tions have also been analyzed for gels with permanent 
bonds [53j or low density gels [54j . 


A. Non-Gaussian parameter 02 (f) 

We begin by looking at how the non-Gaussian param¬ 
eter varies with changing connectivities R, either when 
we go from the liquid to the gel phase, at cj> = 0.50, or 
when we are in a strongly glassy region, <f> = 0.70. For 
the case Tn n k = 10 2 the data is shown in Fig. [T2J 

For <f> = 0.50, Fig. [T2h . a 2 (t) shows at t ~ 1 a peak 
if R = 0 and a shoulder for R > 0. This feature is 
thus related to the dynamical heterogeneity due to the 
steric hindrance mechanism, as in usual glass-forming 
systems [44| . For R > 0 we find in addition a very promi¬ 
nent peak the location of which is basically independent 
of R , which shows that for this packing fraction the time 
at which the system is maximally heterogeneous does not 
depend on R. This observation is in tune with our earlier 
discussion that the timescale for structural relaxation at 
small values of <j> does not change with R if the number 
of bonds is not very large (see Fig. @Ji). The maximum 
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occurs around t ss 10, a time which corresponds to a 
timescale for which a significant number of particles have 
broken the bonds with their neighbors and exit the con¬ 
straints of the network. The location of the peak is thus 
somewhat larger than Tank- As has been documented 
in Ref. at small packing fractions one has on this 
timescale two families of particles, one for the mobile 
particles and the other related to those that are immo¬ 
bile and as a consequence the shape of G s (r, t) deviates 
strongly from a Gaussian. Figure [T2h indicates thus that 
the same behavior persists to the larger volume fraction 
of 0.50. One interesting feature is that the peak height is 
non-monotonous in R in that it increases till R = 6 and 
then decreases again for larger R values. The reason for 
the growth is that an increasing R allows for more diverse 
values of the connectivity for the particles, and hence to 
a stronger variety in the dynamical behavior. On the 
other hand, if R is very large most of the particles are 
strongly connected all the times and hence show a much 
smaller variation in their relaxation dynamics, i.e. they 
behave like in a mean-field like regime. 

Next we study the behavior for cj> = 0.70, Fig. IT2b . 
When there are no bonds, R = 0, we see that «2 (t) is 
peaked at around t = 8. This peak, which is related to 
the steric hindrance, corresponds thus to the one seen for 
( f > = 0.50 at £ « 1 and which, due to the higher density 
has shifted to larger times. For R = 2 (when the perco¬ 
lating cluster of the bonds develops), the location of this 
peak slightly shifts to larger times and one sees the ap¬ 
pearance of a second peak at t ps 2000 (corresponding to 
the timescales for renewal of the network connectivities). 
However, the dominant heterogeneity is still due to the 
local steric hindrances. As R is increased, the location of 
the first peak continues to shift to longer timescales, in 
track with the increasing structural relaxation timescales 
(see Fig. [4J:), and also its height increases, in qualita¬ 
tive agreement with the behavior found in simple glass- 
formers if the coupling is increased [44j |. For R = 6, the 
two peaks have merged, since the timescales for the two 
sources of heterogeneity are nearly the same, and thus 
02 (t) has a single peak. If R is increased even more, the 
position of the peak moves to larger times, but its height 
starts to decrease. The reason for this decrease is likely 
the same as the one we indicated when we discussed the 
data for </> = 0.50, i.e. that for large R the system starts 
to become mean-field like and hence heterogeneities are 
suppressed. 

We also note that the height of the peak is significantly 
smaller than the one for <f> = 0.50. Thus, at large density, 
steric hindrance dominates and even the faster particles 
have less space to move around resulting in significantly 
less non-Gaussian shapes for G s (r, t). As a consequence, 
the dynamics is more homogeneous in this regime com¬ 
pared to the gel at lower cj). 

To elucidate the origin of the two peaks in the non- 
Gaussian parameter at R = 2, <f> = 0.70, we determine 
how the distribution of particle displacements G s (r, t) 
evolves with time. This is shown in Fig. 1131 At short 



FIG. 13: (Color online) Self part of the van Hove function for 
R. = 2 and (j> = 0.70, parameters for which a. 2 (t) shows two 
peaks (see Fig. 112b 1 ). 


times, particle motion is restricted by the local cage of 
neighbors and thus G s {r,t) is a Gaussian at small r. At 
larger distances the distribution has an exponential tail 
which is due to a few particles that have escaped the 
steric hindrance cage, in agreement with the usual het¬ 
erogeneous glassy dynamics in supercooled systems Q. 
The presence of these two processes gives rise to the max¬ 
imum in «2 at short times. At around t = 400, most of 
the particles have escaped from this steric hindrance cage 
and thus the dynamics of the system becomes more ho¬ 
mogeneous with G s (r, t) assuming a Gaussian form and 
hence «2 decreases again. But with time, the motion of 
the particles is again restricted, this time by the bonds 
which have not yet relaxed. Thus again G s (r, t) devel¬ 
ops a Gaussian shape (with a larger width, determined 
by the length of the connecting bonds) and an exponen¬ 
tial tail, implying an increase in 0:2- Later, the network 
eventually relaxes, the particles are diffusive and thus 02 
decreases again. 

The presence of these two different contributions to 
the shape of 02 (t) can be further clarified by varying the 
lifetime of the network, i.e. rii n k, and in Fig.[lT]we show 
this for the case <p = 0.65, R = 6. Since we scan a large 
span of timescales and <22 changes strongly, the data is 
shown in logarithmic scales. For 7ii n k = 1, we see only 
a small peak at t ps 10, which corresponds to the one 
due to the steric hindrance. For Tii n k = 10 2 , we find that 
this peak has increased a bit in height and its location 
has also shifted to t « 30. This change is a consequence 
of the increased effective coupling between the particles 
due to the increased lifetime of the bonds. In addition, 
we see the development of a new (weak) peak at t « 10 3 , 
which is caused by the network renewal process. Now, 
if we increase r^nk further, the location of the first peak 
remains at the same place, since the local crowding ef¬ 
fects are unaffected by the network dynamics. In contrast 
to this, the location of the second peak shifts to longer 
timescales and its height increases strongly. Since the 
position of this peak scales with rii n k, we recognize that 
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FIG. 14: (Color online) Time dependence of a- 2 , for different 
values of riink at R = 6, <j> = 0.65. Inset: For the case of fixed 
bonds (rii n k = oo), a ?2 for the full system (maroon) along 
with those for the clustered (green open symbols) and non¬ 
clustered particles (orange filled symbols). 

this peak is related to the renewal process. Finally, if we 
take the case of permanent bonds between the particles, 
the curve is nearly identical to that for Tn n k = 10 4 , except 
that we see no signature of it eventually decreasing with 
time. The reason for this runaway effects is the fact that 
with increasing lifetime of the network, the unattached 
particles diffuse away and travel long distances before 
the network is again renewed. This results in very long 
extended tails in G s (r,t) which shows up as very large 
values for the non Gaussian parameter. For the case that 
"Hink diverges the dynamics never becomes Gaussian and 
0:2 (t) diverges at long times. 

It is interesting that if one computes, for rii„k = 00, 
0:2 (f) separately for the clustered and non-clustered par¬ 
ticles, only the peak at short times is observed, i.e. the 
one which is due to the crowding effects. This is shown in 
the Inset of Fig. [14] for the case of the system of particles 
where there is a permanently frozen percolating cluster 
and few unattached particles. In this case, the 02 it) for 
the clustered particles follows the curve for the full sys¬ 
tem, shows the bump at short times and becomes at long 
times a constant (see the corresponding data for A 2 (t) 
in Fig. [9]). The result that at long times «2 (t) does not 
go to zero, as it would be the case for vibrations in a 
typical amorphous solid with R = 0, is related to the 
fact that the spanning cluster is a disordered network of 
particles that have different local environments. Thus 
the large variety of slow floppy motions, related to the 
slow relaxation in F s (q,t ) over long times (see Fig. [9]), 
gives a non-Gaussian shape for G s (r, t) even at very long 
times. This non-Gaussianity is also the reason for the ex¬ 
ponential shape of fq N (q), as observed earlier, see Fig. [8] 
For the non-clustered particles, a .2 (i) does show a slight 
decrease beyond the short-time maximum which shows 
that the relaxation dynamics starts to become a bit more 
homogeneous. However, it cannot be expected that the 
non-Gaussian parameter will go to zero even at very long 





FIG. 15: (Color online) Time dependence of \4 for Tu n k = 
100. (a) R = 0 and different values of <j>. (b) <f> = 0.70 and 
different values of R (as marked on plot) (c) R = 4, tj> = 0.65 
for different values of rank (as marked on plot). 

times, since the mobile clusters do have a spread in size 
and hence a different diffusion constant. 

B. Four point susceptibility £) 

We conclude by measuring to what extent the observed 
heterogeneous dynamics is related to a correlated motion 
of the particles. This can be quantified by Xi (<7, i)> de- 
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fined above, measured for q = 6.0, i.e. we look for relax¬ 
ations over distances which are slightly larger than the 
average particle diameter. 

To have a reference we consider first the case R = 0, 
i.e. when there are no polymers and the system is just 
a standard glass-forming system. The time dependence 
of \4 (<?, t) is shown in Fig. 1 1 5b . for different packing frac¬ 
tions. In a gre ement with earlier studies on similar sys¬ 
tems, Ref. [50(, we find that shows a maximum 

at a time t that increases with <j) and which tracks the in¬ 
creasing relaxation time tsh- The fact that the height of 
the peak increases with increasing <f> indicates that the re¬ 
laxation dynamics become more cooperative, also this in 
agreement with previous studies for other glass-forming 
systems [52j . 

Next, we contrast this behavior with the case of in¬ 
creasing gelation, i.e. we fix a volume fraction and in¬ 
crease the number of polymer bonds. For this, we choose 
(j) = 0.70, i.e. the same packing fraction for which we have 
studied the correlation functions and the non-Gaussian 
parameter, since for this </> the effects due to the steric 
hindrance as well as the network constraints are clearly 
observed. The corresponding X4(<Z, t) is shown in Fig. 115b 
for rii n k = 10 2 . For 0 < R < 6, the location of the peak 
in X4(g,t) shifts to larger times and the height of the 
maximum increases. Thus this trend follows the behav¬ 
ior observed in the data for F s (q,t ) as well as 0:2(t), in¬ 
dicating that steric hindrance dominates the relaxation 
process and that the dynamics becomes more and more 
collective. However, for R > 6, the peak height decreases 
again and its location becomes independent of R. We can 
thus conclude that in this regime the relaxation dynamics 
becomes less cooperative, a result that makes sense since 
the rewiring of the network is not really a collective pro¬ 
cess. This result is also in qualitative agreement with the 
decrease of the maximum in 0:2 (f), see Fig. 112b . which 
showed that the dynamics becomes more homogeneous. 

We conclude by investigating how the lifetime of the 
network-bonds influences the function X4(9,t)- For this 
we have calculated this observable for the case <f> = 0.65 
and R = 6, i.e. at intermediate density and in the gel 
phase for which one sees a significant rii n k-dependence in 
the relaxation behavior (see Fig. [7t>). We see, Fig. fl5b . 
that at small and intermediate Ti; n k the location of the 
peak in X4(?,t) moves to larger time and that its height 
increases somewhat, i.e. a behavior that is directly re¬ 
lated to the steric hindrance mechanism in which the ef¬ 
fective cage becomes stiffer due to the increased lifetime 
of the bonds. However, once Tu n k exceeds 10 4 , we see 
that X4(<b^) n0 longer depends on this lifetime, i.e. the 
rewiring of the network no longer affects the cooperativ- 
ity of the dynamics and the latter is solely dependent on 
the steric hindrance. 


VII. SUMMARY AND CONCLUSIONS 

In this paper, we have studied a coarse-grained model 
for a transient network fluid, a system that can be real¬ 
ized in experiments by a oil-in-water emulsion. By tun¬ 
ing the volume fraction (j> of the constituent particles, 
the number of bonds R between the mesoparticle, and 
the lifetime rii n k of these connections, we have scanned 
across the phase diagram to investigate the relaxation 
dynamics of this system, in particular the interplay be¬ 
tween the gel-transition and the glass transition induced 
by steric hindrance effects. 

By analyzing the mean-squared displacement and self 
intermediate scattering function we find that the nature 
of the slowing down depends on the packing fraction: At 
low <f> an increase of R leads to a two step relaxation 
with a plateau that increases continuously with R and 
an a—relaxation time that is independent of R and that 
is related to the lifetime of the bonds. In contrast to this 
we find at large (f) a two step relaxation with a plateau 
height that is basically independent of R whereas the re¬ 
laxation time does depend on the density of bonds. Thus 
these two different behaviors show that the system can 
show a glassy dynamics that is related on one hand to a 
gel transition and on the other hand to a glass transition 
associated with the steric hindrance mechanism. For in¬ 
termediate densities and a certain range of R and ru n k 
the interplay between these two mechanisms leads to a 
decay of the time correlation function that is logarithmic 
over several decades in time, whereas for other combina¬ 
tions it gives rise to a three step relaxation. We note here 
that a similar multistep relaxation pattern has been re¬ 
cently reported for mixtures of multiarm telechelic poly¬ 
mers and oil-in-water microemulsions [55| , suggesting the 
existence of a competition of different arrest length- and 
time-scales in these systems. We have also studied how 
the height of the plateau at long times, i.e. the Debye- 
Waller factor, depends on the wave-vector and found that 
it decays in an exponential manner in q , i.e. very differ¬ 
ent from the Gaussian decay observed in more standard 
glass-forming systems. 

Furthermore we have determined the relaxation times 
of the system and find that these can approximately 
be factorized into a function that depends strongly on 
R and a Vogel-Fulcher type dependence on (j). The 
R —dependent factor shows for R < 10 a strong expo¬ 
nential dependence that is related to the escape of the 
particles from the local cage, whereas for larger values of 
R the dependence is weaker and linked to an Arrhenius- 
like process for bond-breaking. 

By studying the non-Gaussian parameter we have 
probed to what extent the relaxation dynamics of the 
system is heterogeneous. At low packing fraction this 
dynamics becomes extremely heterogeneous, if R is not 
too large, since some of the particles are very strongly 
connected to their neighbors (and hence are immobile) 
whereas others can move almost freely. However, if R 
becomes larger than 6, the relaxation dynamics becomes 
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again quite homogeneous since at any time all particles 
are well connected to their neighbors. At large pack¬ 
ing fractions the non-Gaussian parameter shows a double 
peak structure and, by monitoring the van Hove function, 
we can show that this feature is directly related to the 
two relaxation processes, i.e. the steric hindrance and 
the connectivity of the network. 

Whether or not the relaxation dynamics is cooperative 
can be characterized by the four-point correlation func¬ 
tion Xi{<L t)- We find that for large packing fractions the 
height of the peak in (q,t) increases rapidly with R, 
showing that the strengthened coupling leads to an en¬ 
hanced cooperative motion. However, if R is increased 
beyond 6, this cooperativity decreases again, since the re¬ 
laxation dynamics is strongly dominated by the rewiring 
of the network, i.e. a non-cooperative process. 

Summarizing we see that the interplay between the two 
different mechanisms giving rise to glassy behavior can 
lead a rather complex and unusual relaxation dynam¬ 
ics. In the present study we have focused on a system 
in which the particles are connected by polymers having 
a fixed extension length £ = 3.5 a. For much smaller ex¬ 
tension lengths, one can expect to recover the re-entrant 
scenario observed in colloidal gels Q- For longer exten¬ 
sion lengths, the particles will instead have more space 
to explore and it will likely result in a even more floppy 


network. In the future, it would certainly be worthwhile 
to explore also systems that have polymers with different 
value of t, since this implies multiple localization lengths 
and thus provide novel relaxation scenarios that can sub¬ 
sequently also studied in experiments. And finally we 
recall that our model is motivated by an experimental 
system which has highly nontrivial rheological proper¬ 
ties, e.g. this material flows like a liquid but eventually 
breaks as a brittle solid [56:, (57j, being capable of self- 
healing through thermal fluctuations. Thus, in the future 
we plan to study the rheological properties of our model 
in order to understand the microscopic mechanisms that 
lead to the experimental observations. 
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